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Abstract 

This article presents a derivation of analytical predictions for steady-state dis- 
tributions of netto time gaps among clusters of vehicles moving inside a traffic 
stream. Using the thermodynamic socio-physical traffic model with short-ranged 
repulsion between particles (originally introduced in [Physica A 333 (2004) 370]) 
we firstly derive the time-clearance distribution in the model. Consecutively, the 
statistical distributions for the so-called time multi-clearances are calculated by 
means of theory of functional convolutions. Moreover, all the theoretical sur- 
mises used during the above-mentioned calculations are proved by the statistical 
analysis of traffic data. The mathematical predictions acquired in this paper are 
thoroughly compared with relevant empirical quantities and viewed in the context 
of traffic theory. 
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1. Introduction 

The most of scientific works dealing with traffic theory do not omit to dis- 
cuss certain regularities occurring in traffic micro structure. Explorations of traffic 
micro-quantities (i.e. quantities belonging to an individual vehicle) and efforts 
to predict their statistics are as old as traffic research itself. Virtually, all im- 
portant reviews on traffic science (for example [1, 2, 3], or [4]) try to explain at 
least some basic knowledge on statistical distributions of time-intervals or distance 
gaps among moving cars. Although the theoretical prediction of time-evolution 
for headway-distributions is still extremely vague, some partial achievements have 
been reached in the last years (for example [5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16], 
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or [17]). We try to pick up the threads of those results and get closer to the heart 
of the matter. 

This contribution is focused predominantly on the time-evolution of clear time 
intervals (so-called time-clearances) between succeeding cars passing a given 
point (a traffic detector, typically) located at an expressway. Investigations the 
traffic clearances are significantly advantageous (contrary to explorations of dis- 
tance gaps) because of their direct measurability. Indeed, the most of traffic de- 
tectors gauges the time of vehicle's passage directly, which means that those data 
do not show any systematic or mediated errors (provided the detector is not dam- 
aged). Except the time-clearances we concentrate our attention on the so-called 
multi-clearances, i.e. cumulated time-clearances among n consecutive cars as 
well. Again, such a quantity is (similarly to time-clearance itself) directly measur- 
able, which opens a possibility for detailed statistical analyses of relevant multi- 
clearance evolution (i.e. investigation of multi-clearance distributions with respect 
to the location of traffic systems in the fundamental diagram). 

In fact, the empirical and/or theoretical investigations of traffic multi-headways 
are not sporadic in the physics of traffic. Some of the previous scientific works, 
e.g. [18, 19, 20, 21], or [22], briefly analyze the time multi-clearance distributions 
(predominantly from the empirical point of view) or their statistical variances. The 
reasons for investigations of such a type are obvious, since deeper comprehension 
of changes in traffic microstructure will provide a more thorough insight into the 
convoluted traffic interactions. 



2. Preliminary evaluation of empirical observations 

The vehicular-data records analyzed in this article are typically of the follow- 
ing types. The set 

7f n) = {4 n) G R + | k = 1,2,..., N}, {£ £ {0, 1,2}) 

includes the chronologically-ordered times when the front bumper of kth car has 
intersected the detector- line (located at £th lane of an expressway). The respec- 
tive times r^ ut) when the rear bumper of kth car has intersected the detector-line 
are summarized in the set 7^ out) . Analogously, the set of velocities of individual 
vehicles recorded by the detector is 

V { = [v kc 6 R + | k = 0, 1, 2, . . . , N], (€e {0, 1,2}). 



2 



The associated lengths of vehicles are denoted as du and summarized in the set 
D(. We remark that the indexes i = 0, 1, 2 correspond to the slow, main, and over- 
taking lanes of a freeway, respectively. The slow lane is intended predominantly 
for long vehicles (trucks or lorries) and therefore the respective data will not be 
considered for this research. On contrary, the drivers use the third lane (£ = 2) par- 
ticularly if they are overtaking (and after overtaking-manoeuvre they turn back in 
the main lane) or if they are significantly faster that other cars. All of the above- 
mentioned quantities are (for purposes of this article) considered to be primary, 
which means that they are directly detecable by the traffic detectors. 
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Figure 1: Sample-mean clearance of subsequent vehicles. The average value of real (non-scaled) 
space-clearances among cars is decreasing function of traffic density. The difference between lanes 
is fractional. 

Except the primary traffic quantities we now introduce some important sec- 
ondary quantities, whose values are obtained vicariously, i.e. they are not included 
in the traffic-detector's records. Among the secondary traffic micro-quantities they 
use to be accented the time headways 



r (in) ^(in) 

and time clearances 



Zkt '■- T k{ T (k-\),e 



t . _ _ (out) 

1-kt ■- l kf 



Although both of them are not explicitly included in traffic-data files, their val- 
ues are not burdened with any additional error. Above that, brutto space-gaps 
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between successive vehicles (usually called as the distance headways) are tradi- 
tionally approximated by the relation su '■= Vu Zu that presupposes the constant 
velocity v(r) = Vu during a time period when r e [tSwJ^ 1 ]. As well known, 
such a precondition is questionable, especially in the region of small traffic densi- 
ties where the time headways are too large. However, the influence of a possible 
error is of marginal importance, as apparent from the fact that the headway dis- 
tributions analyzed in small-density regions do not show any noticeable deviation 
from exponential distribution expected for infrequent events (see [11, 2], or [6]). 
Analogously, the distance-clearance is calculated via r# '■= vu hi and represents 
the estimated netto distance between fcth car and its predecessor. Contrary to the 
time headways (clearances) the distance headways (clearances) are burdened by 
the systematic error that is discussed above. Such a error devalues the knowledge 
on distance headway distributions and their evolution. 

Denoting the sampling size by m (in this study there is consistently con- 
sidered m = 50) and number of data- sets by M c (which therefore implies that 
Mf = [Nf/mj), where N e = cardm 1 ? | I is fixed}, one acquires the main data- 
samples 

5 ( m ain) = [(^^^4) 6 7f ) X jf* x y iX Dl | 

k = (j-l)m+l,(j-l)m + 2,...,jm A £= l), (1) 
and the secondary data- samples 

S f ast) = ((r?, rf; l \ v ke , d u ) 6 7f 3 x 7f ut) xV 2 x D 2 \ 

k = (i — \)m + 1, (i - \)m + 2,...,im A € = 2j, 

where ; = 1, 2, . . . , M x and i = 1, 2, . . . , M 2 . For each data-sample 5 J main) (or 5 f ast) 
alternatively) we calculate the local flux 



J (out) _ (in) 

1 jm,t l (j-\)m+lf 

and local average velocity v ; - = m" 1 Zf= (y _i )m+ i v^. The local density gj is then 
estimated (as suggested in [2]) via the fluid-dynamic equation 



Such a expression is understood as one of the approximations suitable for es- 
timation of vehicular density. We add that the incorrectness of the definition 
(2) is caused by the mixing of time and spatial averaging. Besides the macro- 
scopic description of each sample S j we now introduce the mean time-clearance 
(or distance-clearance) by means of definitions 



1 \ 1 1 X 1 

m t—i m. t—L 

respectively. 



m * — ' m 

k=(j-l)m+l k=(J-l)m+l 




Figure 2: Mean scaled velocity (and variance) as a function of traffic density. In the main part of 
the figure there is displayed the average value w which is drawn separately for the fast and main 
lanes. The fluctuations of velocity quantified by the velocity-variance VAR(W) = «~' Z£=i(Wjfc — w) 2 
depending on traffic density are visualized in the inset. 

Since aiming to investigate the essential properties of micro-distributions (i.e. 
statistical distributions of microscopic traffic quantities) we eliminate (in the next 
part of this text) the global trends in those distributions. Predominantly, we elim- 
inate the changes of the average headways caused by the varying traffic den- 
sity. Such an approach corresponds to the techniques examined in the articles 
([8, 11, 9, 21], or [17]), and allows more sophisticated comparison for different 
traffic regimes or for traffic data originated from different countries. The proce- 
dure of the headway re-scaling (see the text below) represents in fact the trivial 
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variant of the so-called Savitzky-Golay smoothing filter (for details please see the 
Ref. [23]) applied to the matrix spectra in the Random Matrix Theory (see [24]), 
for example. Thus, we define the scaled distance-clearances (for the sample S j) 
as 

Xu = — for all k e {(j - l)m + 1,0'- l)m + 2, . . . , jm\ 

r j 

and the scaled velocities as 

w u = — for all ke{(j-l)m+l,(j-l)m + 2,...,jm}. 

r j 

It implies that the mean clearance in each sample is re-scaled to the unit. However, 
for completeness of these deliberations we plot (in the Fig. 1) the dependence of 
the mean distance-clearance on traffic density. Similarly, we analyze the changes 
(averages and standard deviations) of the scaled velocities wu in the Fig. 2. 




Figure 3: Correlation coefficient as a function of traffic density. The value a> wx quantifies a statis- 
tical dependence between the scaled vehicular velocities Wkc and scaled gaps xu to the preceding 
car (in the regions of fixed traffic densities). The squares represent cars moving in the main lane, 
whereas circles correspond to the fast-lane cars. Note that the vehicles in the fast lane show the 
stronger correlations than others. Such a tendency is accented in the regions where the mental 
strain of drivers is stronger (i.e. in the regions of over-saturations). Vertical abscissae demonstrate 
a degree of fluctuations (quantified by the standard deviation) in co wx for main-lane data analyzed 
from different data sources. 

For intentions of analytical calculations executed in one of the next sections 
it is necessary to explore mutual relations between spatial clearances and car ve- 
locities. The standard way how to inspect those relations is to investigate the 
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correlation coefficient. Denoting X = {x x , x 2 , . . . , x n } and W = {w\, w 2 , ■ ■ ■ , w n } 
the sets of scaied distance-ciearances or scaled vefocities respectively, one can 
calculated the Pearson's correlation coefficient using the relation 

COV(X, W) Z'L(x k -lXw k -w) 
" WX VAR(X) • VAR(W) VZLife-D 2 VZZ=i("*-*) 2 ' 
where w = YIk=\ w k/n. Dividing the complete data-sample into sub-samples (where- 
in the traffic density g is the same, approximately) we demonstate (in the Fig. 3) 
that coefficient co wx is localized in the range [0.05; 0.30] and the corresponding 
values are depending on a traffic regime (either free or congested regime) and 
freeway lane. However, in all density-intervals the Pearson's correlation is re- 
duced which means that the hypothesis on a weak dependence between clearance 
and velocity seems to be legitimate. 

3. Mathematical derivation of time-clearance distribution 

As introduced in [8], justified in [14], and elaborated in [11], one of the pos- 
sible ways how to acquire meaningful predictions for distributions of empirical 
distance-clearances is in using the socio-physical traffic model whose elements 
are repulsed by the short-ranged power-law forces and randomized by the socio- 
physical noise, which is supposed to be of thermal-like nature. The influence of 
such a thermodynamical component to the model can be increased/reduced by the 
socio-physical coefficient f3 € [0, oo) reflecting a mental strain under which the 
drivers are during a given traffic situation. Specifically, in the free traffic regimes 
(where the psychological pressure caused by the traffic situation is weak) the pa- 
rameter /3 is almost zero. On contrary, for the congested phase (where interactions 
among drivers are reinforced) the mental strain coefficient /S is large. For the de- 
tailed changes of /S one can inspect the articles [8, 1 1], or [21]. 

Specifically, we consider dimensionless particles moving along a ring whose 
velocities are W\ , w 2 , • ■ • , w N and mutual distances (between subsequent particles) 
are Xi, x 2 , . . . , x N . Introducing the short-ranged repulsive potential 

N 1 

U(x 1 ,x 2 , ...,x N ) = V — (3) 



and socio-physical hamiltonian (see [11]) 

1 N 

f H(w u w 2 ,...,w N ,x 1 ,x 2 ,.,.,x N ) = - 2^(w k - w d f + U(x u x 2 ,...,x N ) (4) 
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one can derive (for details please see [11]) that velocities of such a ensemble (an- 
alyzed in the steady state) are gaussian-distributed, i.e. the associated probability 
density reads 

1 '»'-'';d> 2 

q(w) = — e (5) 

where o 2 is the second statistical moment (variance) and Wd is the optimal velocity 
of drivers. Analogously, in the articles [11, 21] there has been deduced that the 
scaled distance-clearance distribution reads as 

p(r) = A 0(r) e^ lr & D \ (6) 

where 0(x) is the Heaviside's step-function and the normalization factors are 

D*{3+— , (7) 

A' 1 =2 y ^'K l (2^Dfi). (8) 

As verified in [8, 11, 19] the one-parametric distribution-family (6) is in an satis- 
factory agreement with the clearance distribution observed in the real-road data. 
We add that the noise-parameter ft is related to the corresponding traffic density g. 

With the help of formulas (5) and (6) one can mathematically derive the ana- 
lytical form of probability density rj(t) for clear time-intervals between two con- 
secutive vehicles (particles of a model). Since the joint probability density for 
distance and velocity can be (under the condition on weak distance- speed depen- 
dency - see the Fig. 3 again) predicted as g(w, x) = q(w)p(x), the joint probability 
density for time gaps and velocities reads therefore as h(w, t) = wq{w)p{wt). Then 
the time-clearance distribution represents in fact a marginal density 

T](t) = I h(w, t) dw = J wq(w)p(wt) dw. 

After expanding a function f(w) = wp(wt) into the Taylor's series about the opti- 
mal velocity Wd we acquire 

°° 1 d f / 
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where 



fw + e — " , f (/eN). 



dw e diwtY d(wi) 
Hence 

77 a/ Tj WiOf I (w - w d ) ?(v) dw + 
f=J d(wf)< Jr 

where = - Wd/#(w) dw is £-th central statistical moment with three pre- 
rogatived cases /i = 1, Hi = 0, and /i 2 = cr 2 . The latter represents a statistical 
variance (see also (5)). As well known all the odd central statistical moments 
(associated to the Gauss distribution) are zero and even moments comply with 
equalities /j. 2 e = cr 2i (2£ - 1)!!, which leads to the general formula 

*) = wrtw* + g ^ (^(^'w a + ^(^-') . (9) 
Owing to the facts £ w d p{w d t) dt = 1 and 

Jr\o(wO <9(w0 if 1 / 

the TC-distribution (9) is normalized correctly. In a certain sense (as probably 
understandable from the two previous equations) the first factor in (9) forms a 
leading term of TC-distribution, whereas the sum in (9) represents a perturbation 
term only. Above that, the weight of the £th perturbation summands decreases 
with a factor cr 2t /(2£)\\. For practical applications (especially for traffic applica- 
tions) it seems therefore to be meaningful to use an approximate expansion 

17(f) ~ w d p{w d t) + °— -zr^(w d t)t 2 w d + 2—^-(w d t)t) . (10) 
2 \o(wty o(wt) ) 

Such an approximation is legitimate because of extremely low value of velocity 
variance mainly (see the inset in the Fig. 2). Note that the factors cr 2c l(2£)\\ 
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measured in real-road data rapidly vanishes since the maximal value of velocity 
variance is ~ 0.02 in the main lane or « 0.06 in the fast lane. The expansion (10) 
can be re-arranged (after re- scaling clearances) in the formula 

17(f) * p{t) + cr 2 p(t) (fL + !L-- D(t + £)) . (11) 

The influence of the low velocity-variance to TC-distribution is demonstrated in 
the Fig. 4 where the curves (11) are compared to the zero-approximation 

?7(O«A0(Oe^e- Df , (12) 

where the relations (7), (8) hold true. We remind that the normalization and re- 
scaling conditions f n(t) dt = it rj(t) dt = 1 have been applied. 




Scaled Time-Clearance [-] 

Figure 4: Time-clearance distribution (analytically derived). The signs represent the values 
of function (11) calculated for /3 e {0,1/2,3/2,5/2}, whereas the curves display the zero- 
approximation (12) for the same /3. Small discrepancies between signs and corresponding curves 
are caused by fact that statistical variance a 2 of the scaled velocities (measured on freeways) 
is very low. Phenomenologically, this relation can be estimated (for the main-lane data) by the 
inequality a 2 < (45yS/55) 4 exp[-45/3/14]. 

4. Mathematical derivation of multi-clearance distribution 

Now, knowing the one parameter family of the probability densities (11) or 
their zero-approximations (12) one can derive an analytical prediction for the so- 
called nth multi-clearance distribution T n {t) which represents the probability den- 
sity for a clear time gap t among n + 2 neighboring particles. Therefore, the main 
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goal of this section is to quantify (by means of theory of functional convolutions) 
the probability that time-period between two following instants (the first one: the 
back bumper of the kth vehicle is leaving the detector; the second one: the front 
bumper of the (k+n+ l)th vehicle is intersecting the detector line) is ranging in the 
interval [t, t + dt). Using this notation we find that the probability density for the 
time-clearance between two succeeding cars (derived in the previous section) is 
Tj(t) = Uo(t). Regarding the clearances as independent the nth probability density 
r] n (t) can be calculated via recurrent formula 

T] n (t) = IJn-iit) ★ 770(f), 

where symbol ★ represents a convolution of the two probabilities, i.e. 

ln(t) = I T]n-\(s)r]o(t - s)ds. 



As a first resort, we will concentrate our endeavour to studying the low-density 
traffic regimes where the locations of cars are independent and therefore the traf- 
fic flow corresponds (from mathematical point of view) to the Poisson process of 
uncorrelated events. As published in [8, 11, 21] such a situation is characterized 
by the negligible value of socio-physical coefficient /? reflecting a driver's mental 
strain. Thus, the associated clearance distribution reads n(t) = Q(t)e' and sequen- 
tially the multi-clearance distribution for free traffic looks like 

nil r>S2 r*s n -i pS n -\ ,n 

ds„d* l _i...ds 2 dy 1 = 0(O-e" f . 
Jo Jo Jo n\ 

(13) 

As a second resort, we aim to derive a general formula, i.e. formula for non-zero 
J3. For these purposes we use the zero-approximation (12). Applying the method 
of mathematical induction and an approximation of the function 

g n (t,s) = e-^ + ^ *e-^" +l ? (14) 
in the saddle point one can obtain 

rj n (t) = &(t) [ A H - 1 AiT 1 e- f "£e- Ds e-~e- I *- s) ds = 
Jo 

= ®(t)An-iAe Dt f s n - l g n (t, s) ds « 
Jo 

* e(t)A,^Ae^ n+1)2 e- Dt f s"' 1 ds * ®{t)A n -yAn l fz^ l)1 e Dt . 

Jo 
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Hence 

r] n (t) « 0(OA„^e-< (n+1)a e~ D! (15) 
where (after applying the re-normalization procedure) 




n+l 



9C B+1 (2(n+l)V^S). 

This fixes the proper normalization J R rj n (t) dt = 1 . In addition to that the mean 
77th spacing equals to 



f 

Jr 



f ?7„(0 dt = n + \ . 



Note, that (15) holds true also for the limiting case /3 = 0. Indeed, with help of 
formulas lim^o x" +1( K n+ i(x) = (2n)M and rim j8 _ > o + D(/3) = 1 we easily deduce that 




n+l 

9C n+1 (2(n+l)V^)=n!. 

Therefore it holds A n = l/n\, which is in a full consonance with the relation (13). 
Thus, the relation (15) constitutes a zero-approximation for the distribution of 
time multi-clearances. Since we have supposed (in the previous deductions) that 
the variance of scaled velocities is negligible (in a local sense), the form of the dis- 
tribution (15) is a direct consequence of the distribution for spatial clearances and 
the fact, that all cars have practically the same velocity (in a local sense, again). 

Because the empirical measurements show low (but not negligible) variances 
in scaled velocities (see the inset in the Fig. 2), it seems more realistic to derive 
the time multi-clearance distributions under the conditions cr 2 > and cr 2 " « 

for n = 2, 3,4, Thus, we suppose that clear time-intervals among succeeding 

cars are distributed according the rule (11). The detailed analysis of such a rela- 
tions vindicates that the dominating term in the last summand of (1 1) is cr 2 Dtp(t). 
Hence, we surmise that the time clearance follows the low 

r]^\t)«p(t)-cr 2 Dtp(t). 

Here we remark (for mathematical correctness) that this function has to be under- 
stood as an approximative probability density, since the perturbation term cr 2 Dtp(t) 
causes that rf^Xt) does not fulfil the exact mathematical definition. 
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Using the method of mathematical induction we prove below that (under the 
previous surmises) the multi-clearance distribution reads 

rfc\t) * 0<»^Ve~ ( ' l+1)2 < e~ D! (\ - a 2 Dt). (16) 
For completeness, we remark that in the following mathematical operations 

An+2 r>t 

= tt\t) * 1^(0 - 0(0— e~ Dt J e^* Ve"£ ds + 

+ DV 4 0(l)^e- ft f e- ( " +1)2 "a - s)s n+1 e'^ ds - 
n\ Jo 

«! Jo 

-D(r 2 0(O^-e- D? f e- ( ' l+1)2 ^a-5)5"e-- d5 = 
n\ Jo 

= 0(r) ji!!Lrie- ( " +2)2 ^e- Dt (l -o^DO 
(n + 1)! 

there have been used the approximation (14) and the surmise cr 4 « 0. 

5. Empirical multi-clearance distribution vs. analytical prediction 

In this section we will balance the theoretical prognoses deduced in the sec- 
tions 3 and 4 against the empirical multi-clearance distributions of freeway data 
analyzed in the section 2. In the following part of the text we will consider the 
multi-clearances (here: cumulated clearances between five succeeding cars) 

0-l)m+i+3 

*?= Yj f> V=h j = h2,...,M u i=l,2,...,m-3) 

k=(j-l)m+i 1 

enumerated for data sets (1). These multi-clearances are associated with the spe- 
cific traffic density through the formula (2) and re-scaled so that the average 
multi-clearance (in the given data-samples s^ 31 ^) is equal to 4. Now, the multi- 
clearances xf represent a statistical realization of random variable t considered 
in the section 4 and can be therefore confronted with the theoretical probabilities. 
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Figure 5: Statistical distribution of time multi-clearances for low density regimes. The bars shows 
the empirical probability density for clear time-interval among five succeeding vehicles moving in 
main lane (for various density regimes - see legend for details). The blue dashed curves represent 
the prediction (15) plotted for n — 3 and for the fitted value of coefficient /3 = fiat obtained by 
minimizing the weighted error-function (18). Continuous curves (green) display the analytical 
approximation (17) plotted for n — 3 and for the fitted values of /? = ySfi t and s = Sa t specified by 
minimizing the weighted error-function (19). For clearness, we also plot the dash-dotted curves 
(magenta) visualizing the multi-clearance distribution (13) valid for an occurrence of independent 
events/particles. Plus signs illustrates how much the time multi-clearances (gauged by the fast- 
lane-detectors) differ from those detected in the main lane. 
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Figure 6: Statistical distribution of time multi-clearances for high density regimes. The bars shows 
the empirical probability density for clear time-interval among five succeeding vehicles moving 
in the main lane (for various density regimes - see legend for details). The blue dashed curves 
represent the prediction (15) plotted for n — 3 and for the fitted value of coefficient f3 = fiat 
obtained by minimizing the weighted error-function (18). Continuous curves (green) display the 
analytical approximation (17) plotted for n — 3 and for the fitted values of /3 = j8fi t and e = £fi t 
specified by minimizing the weighted error-function (19). For clearness, we also plot the dash- 
dotted curves (magenta) visualizing the multi-clearance distribution (13) valid for an occurrence 
of independent events/particles. Plus signs illustrates how much the time multi-clearances (gauged 
by the fast-lane-detectors) differ from those detected in the main lane. 



15 



To be factual, we divide the entire interval of traffic densities into the subinter- 
vals [0, 3), [1,4), [2, 5), and so on and analyze the freeway multi-clearance distri- 
bution p(x) separately in each subinterval. Such an approach has been examined 
in the articles ([8, 11, 9, 21]) and reflects the known fact that the headway dis- 
tributions are substantially influenced by the changing location of the system in 
the phase diagram. To prevent the mixing of the states with different vigilances 
of car drivers (or with different temperatures of associated heat bath - in the ther- 
modynamic interpretation of vehicular traffic) we make an independent statistical 
analysis of empirical multi-clearances and investigated their evolution with re- 
spect to the traffic density. 

In the Fig. 5 and 6 there are plotted the empirical multi-clearance distributions 
(histograms) against the poissonian distribution (13), zero approximation (15) cal- 
culated for n = 3, and final analytical prediction 

T]i £ \t) = ©<»— rV 16 ^V Dr (l - st). (17) 
6 

The fitted values of the free parameter (J3 in (15), and (3, e in (17)) have been 
determined by means of formulae 

r°°i 2 

/T £it = argmin^o M) \ m (t) - p(t)\ te' t/4 dt, (18) 
Jo 

\rjf(t) - p(t)\ te' t/4 dt, (19) 

o 

i.e. minimizing the statistical distance x(J3, s) = \f(t)-p(t)\ 2 te~ t/4 dt (weighted 
by the factor te~^ 4 ) cumulating the weighted deviations between theoretical pre- 
diction f(t) and empirical frequency p(t). 

As visible in the Fig. 5 and 6 the freeway multi-clearance distribution delin- 
eated for low-density states coincides with the distribution (13), which confirms 
the surmise that vehicles in free traffic regime are moving as independent ele- 
ments. Thus, their statistics is purely poissonian. As the traffic density rises one 
can detect larger deviations from (13), which demonstrates the stronger interac- 
tions among the cars. Roughly speaking, the fitted thermal parameter J3 is increas- 
ing (in both cases: zero-approximation and also the final probability distribution) 
with the traffic density (see the Fig. 7). The regions of the temporal descent in f3 
value agree with the critical regions in the fundamental diagram. To be precise, 
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the transmission between traffic regimes (from free to congested regime) causes 
the transient consolidation of traffic, which leads to a reduction of driver's mental- 
strain. Since (see the [21]) the thermal parameter reflects such a level of mental 
pressure, the detected drop in course of f3 = /3(g) is expectable. We add that the 
evolution of thermal parameter (3 corresponds to the behavior of the associated 
quantity investigated within the scope of the articles [8, 11,21, 26]. Moreover, in 
all traffic states the distribution (17) fits the real-road data more impressively that 
original approximation (15), as comprehensible. Such a fact is clearly visible in 
the Fig. 8 where the statistical distance e) between theoretical and empirical 
distributions are outlined. 




Figure 7: Optimal parameters in theoretical distributions. Circles visualize the value of the fitted 
parameter /?V (see the relation (18)) that minimizes the statistical distance between theoretical 
curve (15) and empirical histogram p(t) (enumerated for main-lane data only). Plus signs corre- 
spond to the value of the fitted parameter jSfit in the relation (19) for which the statistical distance 
Xi/3,s) between the theoretical curve (17) and empirical histogram p{x) is minimal. Similarly, in 
the inset there is plotted the value of the perturbation parameter £fi t considered in (17). 

For illustration purposes (and for better insight into the thermal component 
of the traffic streams) we also plot the dependency of the thermal parameter /?f it 
obtained with help of the formula (19) for separate regions inside the fundamen- 
tal diagram, i.e. we analyze the time multi-clearance distribution in the region 
with a constant traffic flows and constant traffic density. As comprehensible from 
the Fig. 9 the basic trends in (3 = /3(g) dependency fully corresponds to the be- 
havior discussed above. In the region between 35 and 55 veh/km/lane there is 
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transparently demonstrated how the different traffic regimes (i.e. drivers under the 
different mental pressures) influence the statistics of vehicular interactions. 

For completeness, we emphasize again that all the statistical test discussed in 
this article have been executed for the main lane traffic data only. As demonstrated 
in the Fig. 5 and 6 the time multi-clearance distributions detected in the fast-lane 
data show the well-known discrepancies with the main-lane data (predominantly 
in the regions of low densities). For denser streams the fast-lane distributions 
converges to the main-lane distributions, which endorses the belief of scientists 
that in congested traffic states the correlations among the vehicles in different 
traffic lanes are much stronger than in free-flow regime. 
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Figure 8: Statistical distance between theoretical and empirical distributions. We plot the values of 
statistical distances x(fi' £it ) and^(ySfi t , £fit) between main-lane multi-clearance distribution p{t) 
and zero-approximation (15) (pale bars) or final formula (17) (dark bars), respectively. 



6. Summary, conclusion, and future prospects 

To conclude, this article deals with cumulative time clearances among several 
subsequent vehicles passing a given point of an expressway. More specifically, we 
are concentrated predominantly on time intervals between two following occur- 
rences: 1. the rear bumper of Mi car has intersected the detector line, 2. the front 
bumper of (k + n)th car has intersected the detector line. Furthermore, we sub- 
tract the time-periods when vehicles are occupying the detectors, which means 
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that a fundamental quantity for this research is the time multi-clearance. Using 
the local thermodynamical traffic model (originally introduced in [8] and solved 
analytically in [11]) with short-range repulsive potential among the elements we 
have derived (applying theory of functional convolutions) a mathematical formula 
for time multi-clearance distribution. This final theoretical distribution represents 
an one-parametric family of functions, where the one and only free parameter 
reflects the current traffic density. The obtained analytical predictions have been 
successfully compared with statistics of relevant freeway data. The detected cor- 
respondence between theoretical and empirical distributions allows an elaborated 
insight into changes of vehicular-traffic micro structure influenced by the momen- 
tary traffic state. 

As perspicuous from the above-mentioned analysis, the statistical distribution 
of traffic multi-clearances is varying from the Poisson distribution (detected for 
low-density states) to the low-variance distribution (17) that ascertains a pres- 
ence of stronger correlations among neighboring vehicles. The lower statistical 
variance in (17) shows that in the congested traffic states the level of vehicular 
synchronization is much stronger than in free-flow states. Furthermore, the mea- 
sure of such a synchronization can be, as discussed in the article, quantified by the 
thermal-like parameter ft. 




Traffic Density [veh/km/lane] 



Figure 9: Thermal parameter as a function of traffic density and traffic flux. We plot the de- 
pendency of the thermal parameter /3fi t occurring in the formula (19) for various regions of the 
fundamental diagram J = J(g). 
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However, the open problem remains how to approximate the time multi-clear- 
ance distribution for free-flow vehicles moving in fast lanes. Presence of large 
amount of vehicular leaders in free-traffic regimes, i.e. comparable percentage of 
leaders and followers, causes that the clearance distribution is probably a com- 
pound of two partial distributions (first one for leaders, second one for follow- 
ers). Thus, the expected approach leading to analytical predictions for fast-lane 
clearances can be found in the theory of finite mixture distributions (similarly to 
semi-poisson model discussed in [5]). 
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